Microarray experiments: Considerations for experimental design

Microarrays are useful tools to investigate the expression of thousands of genes rapidly. Some researchers remain reluctant to use the technology, however, largely because of its expense. Careful design of a microarray experiment is key to generating cost-effective results. This article explores issues that researchers face when embarking on a microarray experiment for the first time. These include decisions about which microarray platform is available for the organism of interest, the degree of replication (biological and technical) needed and which design (direct or indirect, loop or balanced block) is suitable

Maize microarray annotation database

<p>Abstract</p> <p>Background</p> <p>Microarray technology has matured over the past fifteen years into a cost-effective solution with established data analysis protocols for global gene expression profiling. The Agilent-016047 maize 44 K microarray was custom-designed from EST sequences, but only reporter sequences with EST accession numbers are publicly available. The following information is lacking: (a) reporter - gene model match, (b) number of reporters per gene model, (c) potential for cross hybridization, (d) sense/antisense orientation of reporters, (e) position of reporter on B73 genome sequence (for eQTL studies), and (f) functional annotations of genes represented by reporters. To address this, we developed a strategy to annotate the Agilent-016047 maize microarray, and built a publicly accessible annotation database.</p> <p>Description</p> <p>Genomic annotation of the 42,034 reporters on the Agilent-016047 maize microarray was based on BLASTN results of the 60-mer reporter sequences and their corresponding ESTs against the maize B73 RefGen v2 "Working Gene Set" (WGS) predicted transcripts and the genome sequence. The agreement between the EST, WGS transcript and gDNA BLASTN results were used to assign the reporters into six genomic annotation groups. These annotation groups were: (i) "annotation by sense gene model" (23,668 reporters), (ii) "annotation by antisense gene model" (4,330); (iii) "annotation by gDNA" without a WGS transcript hit (1,549); (iv) "annotation by EST", in which case the EST from which the reporter was designed, but not the reporter itself, has a WGS transcript hit (3,390); (v) "ambiguous annotation" (2,608); and (vi) "inconclusive annotation" (6,489). Functional annotations of reporters were obtained by BLASTX and Blast2GO analysis of corresponding WGS transcripts against GenBank.</p> <p>The annotations are available in the Maize Microarray Annotation Database <url>http://MaizeArrayAnnot.bi.up.ac.za/</url>, as well as through a GBrowse annotation file that can be uploaded to the MaizeGDB genome browser as a custom track.</p> <p>The database was used to re-annotate lists of differentially expressed genes reported in case studies of published work using the Agilent-016047 maize microarray. Up to 85% of reporters in each list could be annotated with confidence by a single gene model, however up to 10% of reporters had ambiguous annotations. Overall, more than 57% of reporters gave a measurable signal in tissues as diverse as anthers and leaves.</p> <p>Conclusions</p> <p>The Maize Microarray Annotation Database will assist users of the Agilent-016047 maize microarray in (i) refining gene lists for global expression analysis, and (ii) confirming the annotation of candidate genes before functional studies.</p

SSHscreen and SSHdb, generic software for microarray based gene discovery: application to the stress response in cowpea

<p>Abstract</p> <p>Background</p> <p>Suppression subtractive hybridization is a popular technique for gene discovery from non-model organisms without an annotated genome sequence, such as cowpea (<it>Vigna unguiculata </it>(L.) Walp). We aimed to use this method to enrich for genes expressed during drought stress in a drought tolerant cowpea line. However, current methods were inefficient in screening libraries and management of the sequence data, and thus there was a need to develop software tools to facilitate the process.</p> <p>Results</p> <p>Forward and reverse cDNA libraries enriched for cowpea drought response genes were screened on microarrays, and the R software package SSHscreen 2.0.1 was developed (i) to normalize the data effectively using spike-in control spot normalization, and (ii) to select clones for sequencing based on the calculation of enrichment ratios with associated statistics. Enrichment ratio 3 values for each clone showed that 62% of the forward library and 34% of the reverse library clones were significantly differentially expressed by drought stress (adjusted p value < 0.05). Enrichment ratio 2 calculations showed that > 88% of the clones in both libraries were derived from rare transcripts in the original tester samples, thus supporting the notion that suppression subtractive hybridization enriches for rare transcripts. A set of 118 clones were chosen for sequencing, and drought-induced cowpea genes were identified, the most interesting encoding a late embryogenesis abundant Lea5 protein, a glutathione S-transferase, a thaumatin, a universal stress protein, and a wound induced protein. A lipid transfer protein and several components of photosynthesis were down-regulated by the drought stress. Reverse transcriptase quantitative PCR confirmed the enrichment ratio values for the selected cowpea genes. SSHdb, a web-accessible database, was developed to manage the clone sequences and combine the SSHscreen data with sequence annotations derived from BLAST and Blast2GO. The self-BLAST function within SSHdb grouped redundant clones together and illustrated that the SSHscreen plots are a useful tool for choosing anonymous clones for sequencing, since redundant clones cluster together on the enrichment ratio plots.</p> <p>Conclusions</p> <p>We developed the SSHscreen-SSHdb software pipeline, which greatly facilitates gene discovery using suppression subtractive hybridization by improving the selection of clones for sequencing after screening the library on a small number of microarrays. Annotation of the sequence information and collaboration was further enhanced through a web-based SSHdb database, and we illustrated this through identification of drought responsive genes from cowpea, which can now be investigated in gene function studies. SSH is a popular and powerful gene discovery tool, and therefore this pipeline will have application for gene discovery in any biological system, particularly non-model organisms. SSHscreen 2.0.1 and a link to SSHdb are available from <url>http://microarray.up.ac.za/SSHscreen</url>.</p

De novo assembly of transcriptomes from a B73 maize line introgressed with a QTL for resistance to gray leaf spot disease reveals a candidate allele of a lectin receptor-like kinase

Gray leaf spot (GLS) disease in maize, caused by the fungus Cercospora zeina, is a threat to maize production globally. Understanding the molecular basis for quantitative resistance to GLS is therefore important for food security. We developed a de novo assembly pipeline to identify candidate maize resistance genes. Near-isogenic maize lines with and without a QTL for GLS resistance on chromosome 10 from inbred CML444 were produced in the inbred B73 background. The B73-QTL line showed a 20% reduction in GLS disease symptoms compared to B73 in the field (p = 0.01). B73-QTL leaf samples from this field experiment conducted under GLS disease pressure were RNA sequenced. The reads that did not map to the B73 or C. zeina genomes were expected to contain novel defense genes and were de novo assembled. A total of 141 protein-coding sequences with B73-like or plant annotations were identified from the B73-QTL plants exposed to C. zeina. To determine whether candidate gene expression was induced by C. zeina, the RNAseq reads from C. zeina-challenged and control leaves were mapped to a master assembly of all of the B73-QTL reads, and differential gene expression analysis was conducted. Combining results from both bioinformatics approaches led to the identification of a likely candidate gene, which was a novel allele of a lectin receptor-like kinase named L-RLK-CML that (i) was induced by C. zeina, (ii) was positioned in the QTL region, and (iii) had functional domains for pathogen perception and defense signal transduction. The 817AA L-RLK-CML protein had 53 amino acid differences from its 818AA counterpart in B73. A second "B73-like" allele of L-RLK was expressed at a low level in B73-QTL. Gene copy-specific RT-qPCR confirmed that the l-rlk-cml transcript was the major product induced four-fold by C. zeina. Several other expressed defense-related candidates were identified, including a wall-associated kinase, two glutathione s-transferases, a chitinase, a glucan beta-glucosidase, a plasmodesmata callose-binding protein, several other receptor-like kinases, and components of calcium signaling, vesicular trafficking, and ethylene biosynthesis. This work presents a bioinformatics protocol for gene discovery from de novo assembled transcriptomes and identifies candidate quantitative resistance genes

Seasonal variation in mycoflora associated with asymptomatic maize grain from small-holder farms in two provinces of South Africa

Seed quality plays an important role in the establishment of healthy crop stands. The aim of this study was to identify the mycoflora associated with maize grain collected over two growing seasons, one experiencing severe drought, from small-holder farms across KwaZulu-Natal (KZN) and the Eastern Cape (EC), two important provinces with maize producing small-holder farmers in South Africa. Asymptomatic maize ears were collected at harvest during two maize growing seasons from farms located in Hlanganani (KZN), Ntabamhlophe (KZN), KwaNxamalala (KZN), Bizana (EC) and Tabankulu (EC). Maize grain was subjected to seed health tests using the agar plate method. The percentage incidence of fungal species isolated from maize grain was determined with species identities confirmed by ITS sequencing. Nine fungal genera were identified with Fusarium species and Stenocarpella maydis the most prevalent. Fusarium verticillioides, Fusarium graminearum and S. maydis were isolated from all sites in both seasons. No fungal species exhibited a higher incidence in the drought season across all sites.  F. graminearum and S. maydis had higher incidences in the wetter season at four and three sites, respectively. F. verticillioides had a greater incidence at the EC sites, particularly the coastal Bizana site. We conclude that local factors have a greater impact than the drought season on the population structure of ear-rot pathogens.  The widespread presence of fungi that are potentially mycotoxin-producing in asymptomatic maize grain poses health risks to consumers and is worthy of further investigation

MADIBA: A web server toolkit for biological interpretation of Plasmodium and plant gene clusters

<p>Abstract</p> <p>Background</p> <p>Microarray technology makes it possible to identify changes in gene expression of an organism, under various conditions. Data mining is thus essential for deducing significant biological information such as the identification of new biological mechanisms or putative drug targets. While many algorithms and software have been developed for analysing gene expression, the extraction of relevant information from experimental data is still a substantial challenge, requiring significant time and skill.</p> <p>Description</p> <p>MADIBA (MicroArray Data Interface for Biological Annotation) facilitates the assignment of biological meaning to gene expression clusters by automating the post-processing stage. A relational database has been designed to store the data from gene to pathway for <it>Plasmodium</it>, rice and <it>Arabidopsis</it>. Tools within the web interface allow rapid analyses for the identification of the Gene Ontology terms relevant to each cluster; visualising the metabolic pathways where the genes are implicated, their genomic localisations, putative common transcriptional regulatory elements in the upstream sequences, and an analysis specific to the organism being studied.</p> <p>Conclusion</p> <p>MADIBA is an integrated, online tool that will assist researchers in interpreting their results and understand the meaning of the co-expression of a cluster of genes. Functionality of MADIBA was validated by analysing a number of gene clusters from several published experiments – expression profiling of the <it>Plasmodium </it>life cycle, and salt stress treatments of <it>Arabidopsis </it>and rice. In most of the cases, the same conclusions found by the authors were quickly and easily obtained after analysing the gene clusters with MADIBA. </p

SND2, a NAC transcription factor gene, regulates genes involved in secondary cell wall development in Arabidopsis fibres and increases fibre cell area in Eucalyptus

<p>Abstract</p> <p>Background</p> <p>NAC domain transcription factors initiate secondary cell wall biosynthesis in <it>Arabidopsis </it>fibres and vessels by activating numerous transcriptional regulators and biosynthetic genes. NAC family member <it>SND2 </it>is an indirect target of a principal regulator of fibre secondary cell wall formation, SND1. A previous study showed that overexpression of <it>SND2 </it>produced a fibre cell-specific increase in secondary cell wall thickness in <it>Arabidopsis </it>stems, and that the protein was able to transactivate the <it>cellulose synthase8 </it>(<it>CesA8</it>) promoter. However, the full repertoire of genes regulated by <it>SND2 </it>is unknown, and the effect of its overexpression on cell wall chemistry remains unexplored.</p> <p>Results</p> <p>We overexpressed <it>SND2 </it>in <it>Arabidopsis </it>and analyzed homozygous lines with regards to stem chemistry, biomass and fibre secondary cell wall thickness. A line showing upregulation of <it>CesA8 </it>was selected for transcriptome-wide gene expression profiling. We found evidence for upregulation of biosynthetic genes associated with cellulose, xylan, mannan and lignin polymerization in this line, in agreement with significant co-expression of these genes with native <it>SND2 </it>transcripts according to public microarray repositories. Only minor alterations in cell wall chemistry were detected. Transcription factor <it>MYB103</it>, in addition to <it>SND1</it>, was upregulated in <it>SND2</it>-overexpressing plants, and we detected upregulation of genes encoding components of a signal transduction machinery recently proposed to initiate secondary cell wall formation. Several homozygous T4 and hemizygous T1 transgenic lines with pronounced <it>SND2 </it>overexpression levels revealed a negative impact on fibre wall deposition, which may be indirectly attributable to excessive overexpression rather than co-suppression. Conversely, overexpression of <it>SND2 </it>in <it>Eucalyptus </it>stems led to increased fibre cross-sectional cell area.</p> <p>Conclusions</p> <p>This study supports a function for <it>SND2 </it>in the regulation of cellulose and hemicellulose biosynthetic genes in addition of those involved in lignin polymerization and signalling. SND2 seems to occupy a subordinate but central tier in the secondary cell wall transcriptional network. Our results reveal phenotypic differences in the effect of <it>SND2 </it>overexpression between woody and herbaceous stems and emphasize the importance of expression thresholds in transcription factor studies.</p

Dual RNA-Sequencing of Eucalyptus nitens during Phytophthora cinnamomi Challenge Reveals Pathogen and Host Factors Influencing Compatibility

Damage caused by Phytophthora cinnamomi Rands remains an important concern on forest tree species. The pathogen causes root and collar rot, stem cankers, and dieback of various economically important Eucalyptus spp. In South Africa, susceptible cold tolerant Eucalyptus plantations have been affected by various Phytophthora spp. with P. cinnamomi considered one of the most virulent. The molecular basis of this compatible interaction is poorly understood. In this study, susceptible Eucalyptus nitens plants were stem inoculated with P. cinnamomi and tissue was harvested five days post inoculation. Dual RNA-sequencing, a technique which allows the concurrent detection of both pathogen and host transcripts during infection, was performed. Approximately 1% of the reads mapped to the draft genome of P. cinnamomi while 78% of the reads mapped to the Eucalyptus grandis genome. The highest expressed P. cinnamomi gene in planta was a putative crinkler effector (CRN1). Phylogenetic analysis indicated the high similarity of this P. cinnamomi CRN1 to that of Phytophthora infestans. Some CRN effectors are known to target host nuclei to suppress defense. In the host, over 1400 genes were significantly differentially expressed in comparison to mock inoculated trees, including suites of pathogenesis related (PR) genes. In particular, a PR-9 peroxidase gene with a high similarity to a Carica papaya PR-9 ortholog previously shown to be suppressed upon infection by Phytophthora palmivora was down-regulated two-fold. This PR-9 gene may represent a cross-species effector target during P. cinnamomi infection. This study identified pathogenicity factors, potential manipulation targets, and attempted host defense mechanisms activated by E. nitens that contributed to the susceptible outcome of the interaction

Search for transient optical counterparts to high-energy IceCube neutrinos with Pan-STARRS1

In order to identify the sources of the observed diffuse high-energy neutrino flux, it is crucial to discover their electromagnetic counterparts. IceCube began releasing alerts for single high-energy ($E > 60$ TeV) neutrino detections with sky localisation regions of order 1 deg radius in 2016. We used Pan-STARRS1 to follow-up five of these alerts during 2016-2017 to search for any optical transients that may be related to the neutrinos. Typically 10-20 faint ($m < 22.5$ mag) extragalactic transients are found within the Pan-STARRS1 footprints and are generally consistent with being unrelated field supernovae (SNe) and AGN. We looked for unusual properties of the detected transients, such as temporal coincidence of explosion epoch with the IceCube timestamp. We found only one transient that had properties worthy of a specific follow-up. In the Pan-STARRS1 imaging for IceCube-160427A (probability to be of astrophysical origin of $\sim$50 %), we found a SN PS16cgx, located at 10.0' from the nominal IceCube direction. Spectroscopic observations of PS16cgx showed that it was an H-poor SN at z = 0.2895. The spectra and light curve resemble some high-energy Type Ic SNe, raising the possibility of a jet driven SN with an explosion epoch temporally coincident with the neutrino detection. However, distinguishing Type Ia and Type Ic SNe at this redshift is notoriously difficult. Based on all available data we conclude that the transient is more likely to be a Type Ia with relatively weak SiII absorption and a fairly normal rest-frame r-band light curve. If, as predicted, there is no high-energy neutrino emission from Type Ia SNe, then PS16cgx must be a random coincidence, and unrelated to the IceCube-160427A. We find no other plausible optical transient for any of the five IceCube events observed down to a 5$\sigma$ limiting magnitude of $m \sim 22$ mag, between 1 day and 25 days after detection.Comment: 20 pages, 6 figures, accepted to A&

Neutrinos below 100 TeV from the southern sky employing refined veto techniques to IceCube data

Many Galactic sources of gamma rays, such as supernova remnants, are expected to produce neutrinos with a typical energy cutoff well below 100 TeV. For the IceCube Neutrino Observatory located at the South Pole, the southern sky, containing the inner part of the Galactic plane and the Galactic Center, is a particularly challenging region at these energies, because of the large background of atmospheric muons. In this paper, we present recent advancements in data selection strategies for track-like muon neutrino events with energies below 100 TeV from the southern sky. The strategies utilize the outer detector regions as veto and features of the signal pattern to reduce the background of atmospheric muons to a level which, for the first time, allows IceCube searching for point-like sources of neutrinos in the southern sky at energies between 100 GeV and several TeV in the muon neutrino charged current channel. No significant clustering of neutrinos above background expectation was observed in four years of data recorded with the completed IceCube detector. Upper limits on the neutrino flux for a number of spectral hypotheses are reported for a list of astrophysical objects in the southern hemisphere.Comment: 19 pages, 17 figures, 2 table